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Abstract. This paper concerns the regularity and geometry of the free boundary in the optimal 
partial transport problem for general cost functions. More specifically, we prove that a locally 
Lipschitz cost implies a rectifiable free boundary and a C 1 cost implies a locally Lipschitz free 
boundary. As an application, we address a problem discussed by Caffarelli and McCann [T] 
regarding cost functions satisfying the Ma-Trudinger-Wang condition (A3) : if the non- negative 
source density is in some L P (R") space for p 6 (2±i,co] and the positive target density is 
bounded away from zero, then we show that the free boundary is a semiconvex C,j" hypersurface 
with a = 2p(2n-i)-n+i ■ Furthermore, a corresponding regularity theory is initiated in the 
Riemannian setting. 



In the optimal partial transport problem, one is given two non-negative functions / = 
fxn, 9 = 9XA £ L (W 1 ) and a number < m < min{||/||£i, H^Hli}. The objective is to 
find an optimal transference plan between / and g with mass m. A transference plan refers to a 
non- negative, finite Borel measure 7 on W 1 xl", whose first and second marginals are controlled 
by / and g respectively: for any Borel set A C M n , 



where c is a non-negative cost function. 

Issues of existence, uniqueness, and regularity of optimal transference plans have recently been 
addressed by Caffarelli & McCann PQ, Figalli [2], [3], and Indrei [5]. Indeed, existence follows 
readily by standard methods in the calculus of variations. However, in general, minimizers fail 
to be unique and it is not difficult to construct examples when \supp(f Ag)\ > (with | • | being 
the Lebesgue measure and supp(f A g) the support of / Ag := min{/, g}). Nevertheless, Figalli 
proved that under suitable assumptions on the cost function, minimizers are unique for 



Proposition 2.2 and Theorem 2.10]. Up to now, the regularity theory has only been developed 
for the quadratic cost. In this case, if the domains f2 and A are bounded, strictly convex, and 
separated by a hyperplane, Caffarelli and McCann proved (under suitable conditions on the 



1. Introduction 




An optimal transference plan is a minimizer of the functional 




11/ Afli|ii(]R«) < m < min{||/|| L i (R n), ||5|| L i (R n)} 
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initial data) that the free boundaries dU m n Q and dV m n A are locally C 1,a hypersurfaces up to 
a closed singular set S contained at the intersection of free with fixed boundary, see [U Corollary 
7.15]. In the case when there is overlap, Figalli proved that away from the common region finA, 
the free boundaries are locally C 1 [2 Theorem 4.11]; Indrei improved this result by obtaining 
local C > a regularity away from the common region and up to a relatively closed singular set S, 
necessarily contained at the intersection of fixed with free boundary, see [U Corollary 3.13] for 
a precise statement. Moreover, under an additional C ' regularity assumption on Q and A, he 
proved that S is T-L n ~ 2 a- finite and in the disjoint case S C S with T-L n ~ 2 (S) < oo pfl Theorem 
4.9]. 

All of the aforementioned regularity results were developed for the quadratic cost. Our main 
aim in this paper is to obtain free boundary regularity for a general class of cost functions J-q 
satisfying the Ma-Trudinger-Wang (A3) condition introduced in [8J and used in the development 
of a general regularity theory for the potential arising in the optimal transportation problem 



(see Definition 2.4). With this in mind, we prove the following theorem which solves a problem 



discussed by Caffarelli and McCann [TJ Introduction]: 

Theorem 1.1. Let f = fxn £ L p (M. n ) be a non-negative function with p G (^-^,oo] ; and 
9 = 9XA a, positive function bounded away from zero. Moreover, assume that £1 and A are 
bounded, A is relatively c-convex with respect to 0, U A, and 17 n A = 0. Let c G J~o and 
m G (0, min{||/||^i, H^Hli}] .Then dU m H O is locally a C 1,a graph, where dU m HQ is the free 
boundary arising in the optimal partial transport problem and a = 2p(2n~i)-n+i ' 

The strategy of our proof is as follows: first, we utilize a result of Caffarelli and McCann [I] 
to prove that the active region is generated by level sets of the cost function. Thus, the free 
boundary is locally a suprema of these level sets (at least at the heuristic level). Then, thanks to 
the assumptions on c, we prove that the free boundary enjoys a uniform interior cone condition; 
this implies that it is locally a Lipschitz graph in some system of coordinates. To obtain the C, 
result, we connect the free normal with the solution of a generalized Monge- Ampere equation 
for (A3) cost functions, and employ regularity results established by Loeper [6] and refined by 
Liu 0. 

In fact, thanks to the method developed by Figalli [2], o ne ca n localize the problem and 



eliminate the disjointness assumption Q, n A = 0, see Corollary 3.7 As previously discussed, we 



achieve the proof of Theorem 1.1 by first establishing the free boundary to be locally Lipschitz 



(see Proposition 3.2). For this, we only need the cost to be C . We also show that with merely a 



locally Lipschitz assumption, the free boundary is rectifiable (see Proposition 3.1). These weak 
results enable us to also study the partial transport problem on Riemannian manifolds with cost 
equal to the square of the Riemannian distance. 

The rest of the paper is organized as follows: in §2, we state and prove some preliminary 
facts. Then in §3, we proceed with the proof of Theorem [LTJ and in §4 we address the problem 
in a Riemannian setting. 

2. Preliminaries 

Definition 2.1. Given an (to — \) -plane ir in M m , we denote a general cone with respect to it 
by 

C a (vr) := {z G M m : a\P v (z)\ < P^z)}, 
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where it © it 1 - = W" 1 , a > 0, and P 7r (z) & P n ±(z) are the orthogonal projections of z G M. m onto 
it and it 1 - , respectively. 

Definition 2.2. A domain D is said to satisfy the uniform interior cone condition if there exists 
a > and 5 > such that for all x G dD there exists v x G S n_1 so that 

(y + C a (^)) n B s (x) CDD B s (x), 

for all y G D n B§(x). We define the profile of such domains to be the ordered pair (5, a). 

Definition 2.3. A domain D C M. n is said to satisfy a uniform interior ball condition if there 
exists r > such that for all x G dD, there exists v x G S n_1 for which B r (x + ris x ) C D. 

Definition 2.4. We denote by J- , the collection of cost functions c : M. n x IR" — > M that satisfy 
the following three conditions: 

1. c G C 2 (R n x W 1 ); 

2. c(x, y) > anc? c(x, y) = only for x = y; 

3. (Al) Forx,p G M n , there exists a unique y = y(x,p) G M™ sitc/i that X? x c(x , y) = p (left twist); 
similarly, for any y,q G M n , t/iere exists a unique x = x(y, q) G M n smc/i that V y c(x, y) = q 
(right twist). 

Furthermore, we denote by J-q, the set of C 4 (M n x W 1 ) cost functions in T that satisfy: 

4. (A2) det(V (a . )2/ )c) / for all x,y G M n ; 

5. (A3) Forx,p G M n , 

Aij, k i( x ,P)Z£jVkVi > co\t?H 2 V (,i)er,(^) = 0,co>0, 
where A^^i := c r ' k c s,l (c m,n Cij :rn c n:rs — Cij^ rs ), and (c l,jf ) is i/ie inverse matrix of (cij). 
Remark 2.5. Some authors use the notation (A3) s in place of (A3) in condition 5 of Definition 



Definition 2.6. A set V C M n is c-convex with respect to another set U C W 1 if the image 
c x (x, V) is convex for each x G U. 

Lemma 2.7. Let Cl C M. n , Act™ be two domains and c G C 1 (IR n x R n ); assume 

(2.1) b := inf c(x,y) > 0. 

xefi,j/eA 

(2.2) &! := inf |V x c(x,y)| > 0, 

xer2,yeA 

T/ien for any b > bo and y G W 1 , the domain Ey := {x : c(x,y) < b} satisfies a uniform 
interior cone condition with profile depending only on bo, ||c||ci, md the modulus of continuity 
of c x . 
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Proof. Fix y £ A, and consider (ft(x) := c(x,y). Then for a fixed point xo £ {x : (j)(x) = b}, we 
choose a coordinate system such that x n is the direction of the normal to the level set pointing 
into the sublevel set {x : 4>(x) < b} and xq is the origin. Let < 9 < | and note that if x has 
angle 9 with e n , then 

(f)(x) = 0(0) + V0(x) • x + o(x) < 0(0) - 611x1005(6*) + o(x). 

Now since c E C 1 (IR n x W 1 ), by the uniform continuity of c x we have o(x) < ^b\\x\cos(9), for 
x £ -6,5(0) and 5 > (depending on b±, 9, and the modulus of continuity of c x ). Thus, 0(x) < b 
when x has angle at most 9 from e n and is in the 5-ball around the origin. □ 



Remark 2.8. By the positivity ofb\ in (2.2), it follows that we may take v x := — j^^yy as the 



direction of the cone at each point x £ dE b n SI and y £ A. 



y 

Remark 2.9. Note that since c £ C 1 , for a sufficiently small 5 > 0, we may take 9 arbitrarily 
close to I in the proof of Lemma 2.1 . In other words, given any a > 0, there exists 5(a) > 
such that (5(a), a) can be taken as a profile for the level sets of c under the assumptions of 
Lemma 12? 



If the domains O and A have disjoint closures and c(x, y) is a continuous cost function vanishing 
only on the diagonal x = y, then (2.1 ) follows readily. The next lemma gives sufficient conditions 
for ([272]) to hold. 

Lemma 2.10. Let 9, C W 1 , A C M n be two domains and c £ ^(W 1 x R n ); assume 

inf c(x, y) > 0, 

and suppose c satisfies the left twist condition and condition 2 in Definition 2.4. Then 

inf \V x c(x,y)\ > 0. 

Proof. Suppose on the contrary that there exists (x, y) £ £2 x A for which V x c(x, y) = 0. 
Let 4>(x) := c(x,x); using condition 2, 4>(x) > and 4>(x) = only for x = x. Therefore, 
V x c(x, x) = 0, but by uniqueness, we must have x = y (using the left twist condition), and this 
contradicts the positivity of 60 • D 

Lemma 2.11. Let c £ T , and consider two domains O C W 1 , A C K n with disjoint closures; 
set 

Then for any b> bo and y £ M. n , the domain Ey := {x : c(x, y) < b} satisfies a uniform interior 
ball condition with radius r = r(bo,bi, 1 1 c| | ^2 ) > 0, where b\ is defined by (2.2). 

Proof. First, note that since c £ J- we have 61 > by Lemma |2.10 Now for a fixed yo £ A, 
denote 4>(x) := c(x,yo). Then for a fixed point xq £ {x : (j)(x) = b}, we choose a coordinate 
system such that x n is the direction of the normal to the level set pointing into the sublevel set 
{x : 4>(x) < b} and xo is the origin. Now let r := where C2 = 1 1 c| | ^72 , and consider the ball 
B r centered at (0, . . . , r) with radius r. In particular dB r touches the origin. Now we will show 
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that B r C {x : 4>{x) < b}: indeed, it is simple to see that for x G B r , cos(0) > i^J = ^jjp, where 
9 is the angle between x and e n . Therefore, 

tf>[x) < 0(0) +V0(O) >x+ °^\xf 
= b-\V(f>(0)\e n -x + ^\x\ 2 

□ 



Remark 2.12. 5y interchanging the roles of x and y in Lemma 2.11, a similar statement holds 
for E\ := {y : c(x,y) < b}. 

3. Regularity theory 

In this section, we prove various regularity results for the free boundary under minimal as- 
sumptions on the cost function. 

Proposition 3.1. (Rectifiability) Let f = fxn, 9 = 9XA be non-negative integrable functions 
and m G (0,min{||/||^i, H^Hli}] • If c : R n x IR n — > R + is locally Lipschitz in the x variable, and 

(3.1) i 8 x c, 

where d x c is the Clarke subdifferential of c, then the free boundary arising in the optimal partial 
transport problem is (n — l)-rectifiable. 

Proof. First, by utilizing [1, Corollary 2.4] we obtain 

U m f)Q:= [J {x G Q : c(x,y) < c(x,y)}. 

(S,y)S7m 



Next, since c is locally Lipschitz and (3.1) holds, we may apply the nonsmooth implicit function 



theorem [101 Theorem 10.50] to deduce that for a > and y G A, the level set 

{x G R n : c(x, y) = a] 

is locally an (n — l)-dimensional Lipschitz graph. Now, for x G dU m D f2, it follows that 

x G d{x G ft : c(x,y) < c(x,y)}, 
for some (x,y) G 7 m . Hence, there exists a profile (5 x ,a x ) such that 
(3.2) (x + C a >^)) n B 5x (x) c (u m nfi)n B Sx (x), 

for some v x G S n . Consider the sets 

A] -.= [z g (du m n ft) n B Sx (x) :6 e >h,a K < j}, 

and note that by the argument leading to ( |3.2| ), each z G dU m n ft has a profile (6 Z , a z ). Now 
for each j G N, we may select ej > so that P := {yi\iZ{ is a sufficiently fine partition of S n_1 
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in the following sense: for each v G §> n_1 , there exists z/« G P so that \v — vi\ < ej, where ej is 
chosen so that \v — Vi\ < ej implies v G C2j(i / j L ) (note: ej \ as j — > oo). Next, let 

z G Afj := {z£A*: \v z - v*\ < e^}. 

By using v z G C2j(v^~) and a 2 < j, it follows that there exists ay > so that 

Moreover, since <5 Z > j, we may select 5j > and combine it with (3.2) (with x replaced by z) 
to deduce 

(z + c aj {ul)) n b &j (x) c (c/ m n n) n (x). 

Thanks to this cone property, it is not difficult to show that for each i , j G N, Afj is contained on 
the graph of a Lipschitz function (generated by suprema of the cones with fixed opening given 
by ctj). This shows that A*- is (n — l)-rectifiable. Moreover, 

oo 

(du m nn)nB 5x (x) = \J{jAf j . 

j=l i=l 

Next, let (dU m n £l) s '■= {x G dU m H 12 : dist(x, dil) > s}. By compactness, there exists 
{^fclfcil C (dU m n 0) s C 9f/ m n n so that 

n(s) 

(du m nn) s =\J (du m n n). n B 5xk (x k ). 

k=l 

From what we proved, it follows that 

n(s) oo "S- 

(du m nn) s = |J (J 

fc=l j=l i=l 

where each yl?* is (n — l)-rectifiable. Thus, by taking s — > 0, we obtain the result. □ 

Proposition 3.2. (Lipschitz regularity) Let f = fxn, 9 = 9XA be non-negative integrable 
functions and m G (0, min{||/| j^i , | \g\ \l 1 V\ ■ A ssume that A is bounded and c-convex with respect 
to Q. If c G C 1 (M n x W 1 ) and satisfies (2.1) and (2.2), then the free boundary in the optimal 
partial transport problem is locally a Lipschitz graph. 

Proof. By our assumptions, we have 

(3.3) U m nn:= \J {x G n : c(x,y) < c(x,y)}, 

see [U Corollary 2.4]. Next, let x G dU m n 0, and note that since 

x G d{x G ft : c(x,y) < c(x,y)}, 



Lemma 2.7 implies the existence of a profile (5, a) so that 

(x + C a {i^)) n b s (x) c (C/ m n O) n B^x), 
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C x (x,T m {x)) 

\c x (x,T m (x))\ 



[see Remark 



where v x : 

7 m is supported. Note that T m may be multi-valued (i.e. there may exist [x, y\ 



2.8) and T m is the map in [TJ Lemma 2.3] on whose graph 

and (x, 2/2) hi 

the support of 7 m with y\ 7^ y-i). In this case, we select one of them as the value of T m {x) (the 
key fact is: {T m (x)} C A). Now for z G dU m n Q n B$(x), consider the bounded, convex set 
Ca,(z, A) (the boundedness follows from the C 1 regularity of c and the boundedness of Q and the 
convexity is a result of the c-convexity assumption of A) . As 



< inf c x [ 



z,y) 



n-l 



(this follows from (2.2)), the origin is not in the closure of c x (x, A). Thus, we may find £ x € S 
and u x > so that — c x (x, A) C C Ux (^) (i.e. we may trap a bounded, convex set whose closure 
does not containing the origin inside a cone of opening smaller than it). Up to possibly decreasing 
co x slightly, we may assume that the boundary of — c x (x,A) is disjoint from the boundary of 
Cu x (£x)'i hence, by C 1 regularity of c, there exists S x > so that —c x (z,A) C Cu x {£,x) for all 
z G B$ x (x); recall that for z G dU m n B$ x (x), 



(z + c a {v^)) n b s (z) c (c/ m n 0) n £ 5 (z), 

X£x) ( since ^2 = ~ |c*(z'r^(z))| and _c ^( z ' A ) c ^(£3) ); now the angle between 
v z and v is strictly less than 5 due to the fact that the opening of C^^j" ) is strictly less than 



(3.4) 

with v z G 



it. Thus, thanks to Remark 2.9, we may assume £ x G C a (v z ) (by picking a > sufficiently 



small). Hence, there exists a x > a so that (see Figure [T]) 

zeaC/ m nB a;c (a;) 

Combining this information with ( |3.4[ ) yields that for all z G <9t/ m Pi B$ x (x), 

(3.5) (* + c^)) n b,„(x) c (u m n ft) n b 5x (x) 

(we may assume B$ x (x) C B$(z) by choo sing 5 X sufficiently small relative to 5). In fact, by 
possibly taking a x larger and 5 X smaller, (3.5) hol ds fo r all z G U m H Bg x (x): indeed, let z G 
U m n Bs x (x). Then (z,T m (z)) G 7 m and thanks to (3.3), z G ^ with 6 = c(z,T m (z)) > 0. 



Remark 



2.{ 



implies the existence of v z 



^ G S n_1 so that 

(y + c a (^)) n b 5 (z) c E b Tm{z) nncu m nn, 



\c x (z,T m (z))\ 
?b 



for all y G -E^ / a n Bg(x). In particular, 

( 2 + c a (iA))nB^)c[/ m nn. 

Thus, by possibly taking 5 X smaller, if necessary, we may assume Bs x (x) C B$(z), an d if z is 
close enough to x we also have v z G Cu^Ca; )> thus, repeating the argument above from (3.4) to 
(3.5) yields the result. Therefore, we proved the existence of 5 X > 0, a x > 0, and £ x G S n_i so 
that for all z G U m D Bs x (x), 

(3.6) 
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-c x (x,A) 








\ X ^ 








\ \ 




\ \ 







Figure 1. C a M) C r\ z edu m nB 5x ( x ) C a {^)- 



By rotating and translating the coordinate system, we may assume x = 0, £ x = — e n , and 
7r := = M n_1 ; moreover, note that the cone C Qo (7r) is symmetric with respect to the e n axis. 
Define <p : B vo (0) C R n_1 ^ R by 

4>{z') : = sup Ky(z'), 

J/:=(2/',»n)e9C7 m nB„ (0) 

where -B^ o (0) := Proj n (B r)0 (0)) and -Ky is the cone function at the point y on the free boundary 
generated by C ao (ir). Note that 4> is Lipschitz since it is the supremum of Lipschitz functions with 
bounded Lipschitz constant (depending on the opening of the cones). Moreover, by construction 
we have 

(3.7) dU m n B VQ (0) C graph ^(py 

Now we claim that there exist constants d, d G (0, 1) with d depending on the profile of the level 
sets of the cost function, so that 

(3.8) graph^l^o) C dU m n B^Jfij. 

Indeed, pick any d G (0, 1); we may select a constant d = d(d,ao) > small enough, so that 
the graph of 0(2?^ (0)) is contained in 5^(0) (this is possible, since 4> has a uniform Lipschitz 
constant in B Vo (0) which depends only on the profile of the level sets) . Let y G graph 4>\^ , Q) C 

Bi (0). If y £ dU m fl Br (0), then since y is on an open cone with opening inward to U m D $7, 
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it follows that y G U m n fl. Since dU m H -Sj^q) is compact, for > small, it follows that 
Qe(y) H 9J7 m fl Bt = 0, where = Qe{y) is a small cylinder whose interior is centered at y 
and whose base diameter and height is equal to 9; in particular, Qe H graph </>|^ d ^ does not 
contain any free boundary points. Next we consider a general fact: let w G graph <j>\ ^ ^ \ <9?7 m , 
Lt(w) := w + te n , and 

s(w>) := sup t; 

{t>0:L t (w)eu m nn} 

note that since w G graph 

(3.9) s(iy) > s(it)) := sup i, 

{t>0:Lt(tD)£B, (0)} 

(otherwise it would contradict the definition of <f) as a suprema of cones in i?^ (0) and w as a 
point on the graph of <p). Next, keeping the base fixed, we enlarge the height of the cylinder along 
the {y + te n :iGl} axis in a symmetric way (with respect to the plane y n + ir = R n_1 ) so that 
it surpasses At]q; we denote the resulting cylinder by Qg. By (3.9) we have Qg n B vo C U m n fl. 
Then we increase its base diameter, 9, until the first time when Qg hits the free boundary 
dU m n fl inside B„ (0) , and denote the time of first contact by 9 and a resulting point of contact 
by yg (note that since 6 dU m n B Vo (0i), this quantity is well defined). Since 4> is a continuous 
graph in B m (0), and both y and ye are on the graph, we may select a sequence of points 
yk G graph ^ n Qe such that y^ — > ye (by connectedness of graph (t>\§ n ( ) ^ Since 

ye G B vo (0) is an interior point, for /c sufficiently large we will have yk G B vo (0) n Qe, see Figure 
[2] Thus, by definitio n of 9, we will have that the yk are not free boundary points but on the 
graph of cp; thus, by s(y k ) > s(y k ), and this implies y k := y k + s(y k )e n G 8B m (0) C\U m . 

By ( |3.6[ ) we have 

(v* + C Q0 (7r)) n 5^(0) c[/ m nfi. 

However, for large /c, ye G (y k + C ao (n)) (see Figure [2j and this contradicts that yg is a free 



boundary point, thereby establishing (3.8). Thus, combining (3.7) and (3.8) we obtain that in 



a neighborhood around the origin, the free boundary is the graph of the Lipschitz function </>; 
hence, the normal to the graph exists for 'H n ~ l a.e. z' G B Vo (0) and has the representation 
(D<f>{z ), 1)^ ^ point (z',0(z')) where it exists. 

□ 

Corollary 3.3. (Semiconvexity) Let f = fxn and g = gx\ be a non-negative integrable func- 
tions. Assume fl n A = and that A is bounded and c-convex with respect to fl. If c G J 7 , then 
dU m n fl is locally semiconvex. 



Proof. By Proposition 3.2, it follows that dU m nfl is locally a Lipschitz graph, and Lemma [2.11 



implies a uniform interior ball condition with v r = — , Cx { x ''^ n { x }}, as the direction of the ball (c.f. 
Remark |2.8[ ). Thus, locally, the free boundary may be represented as a suprema of spherical 
caps (see [TJ Section 5]), and this readily implies semiconvexity. □ 

Remark 3.4. By a localization argument, one may remove the disjointness assumption in Corol- 
lary 3.3. Indeed, the precise statement (and proof) is similar to Corollary 3.1 below. 

Now we are ready to prove our main result: 
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Figure 2. y e e U m n fJ. 



Proof of Theorem 11.11 First, note that by [21 Remark 2.11], there exists a unique solution to 
the optimal partial transport problem, and it is denoted by the transference plan j m . Moreover, 
by the results of [21 Section 2], this solution has the form 

7 m := (Id x T m ) # f m = (T' 1 x Id) # g m , 

for some invertible map T m , where f m and g m are the first and second marginals of 7 m , respec- 
tively. Now T m is constructed by solving the classical optimal transport problem between the 
densities f + (g — g m ), g + (/ — f m )- Indeed, if we denote this solution by 7, then since c £ J^, 
by applying the classical theory, we know 7 is supported on the graph of a function, say T m . 
Moreover, by [21 Proposition 2.4] and [2, Remark 2.5], it follows that 

(3.10) 7 = Tm + (Id x Id) # ((f - f m ) + (g- g m )), 
and there exists a potential function fy m which satisfies: 

(3.11) V x c(x,T m (x)) = W m (x). 

Now by [21 Theorem 2.6 and Remark 2.11], we have (T m )^(/ m + (5 — <7m)) = 5 (he. T m will not 
move the points in the inactive region). Let /' := f m + (g — g m ) and note that thanks to our 
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assumptions on / and g, 

\det(D 2 xy c)\^— g Lnu m nn). 
y g(T m ) 

Thus, we may apply [5] Theorem 1] to obtain 

(3.12) * m g c^iUm n n), 

Now thanks to Theorem 3.2, we know that the free boundary is lo cally a Lipschitz graph, with 

c x (z,T m (z)) 
\c x {z,T m {z))\' 



T-L n 1 a.e. defined normal v z = — r^f4w^\i • Thus by combining (3.11) and (3.12), we readily 



obtain the result. 

□ 

Remark 3.5. (E xcha nge symmetry) By reverse symmetry, we may interchange the roles of f 



and g in Theorem 1.1 in order to obtain C, regularity of dV m n A. 



Remark 3.6. ( Geometry of c- convex domains) For a geometric description of c-convex domains, 
see [HI Section 2.1-2.3]. For example, based on a calculation therein, one can prove the following: 
suppose A is a bounded, open convex set with smooth boundary and Q C W 1 . Let 

ai := inf I detc xv (x,y)\ > 
igHua.j/gA 

and a2 '■= ||c(-, ■ ) 1 1 c 3 / f or a fixed i g U A, if the principal curvatures of dA are greater than 
then c x {x,A) is convex. 

Corollary 3.7. (N on- disjoint case) Let f = fxn G L p (W l ) be a non-negative function with 
p G (^-^,00], and g = g\A a positive function bounded away from zero. Moreover, assume that 
f2 and A are bounded and A is relatively c-convex with respect to a neighborhood o/OUA. Let 
c G J-q and m G (||/ A g||£i,mm{||/||£i, ||<7||i,i}] • Then away from d(£l n A), it follows that 
dU m n is locally a C 1,a graph, where dU m n Q is the free boundary arising in the optimal 
partial transport problem and a = 2 P (2n-i)-n+i ' 

Proof. Note that by [21 Remark 3.2] and [21 Remark 3.3], we have Q, n A C U m n 0. Therefore, 
the free boundary does not enter the common region f2 n A. Now let x G (dU m nfi)\ d(Q n A). 
Choose r x > so that B Tx (x) n (U m n O) does not intersect A. Thus, dist(B rx (x) n (U m n 
n),T m (B rx (x) n (u m n n)) > o, and so we may apply Lemma 



2.11 to obtain that all level sets 



of the cost function c G To intersecting, say, Bm (x) n (U m n J 2) have a uniform interior ball 
condition (in fact, a uniform interior cone condition is sufficient). Since all of Figalli's results 
used in the proof of Theorem |1.1| are also valid in the non-disjoint case, by localizing the problem 
in this way, we may proceed as in the proof of Theorem |1.1| to obtain the result. □ 



Remark 3.8. The analog of Remark 3.5 is also valid for the case Q n A ^ 



4. Extensions to Riemannian manifolds 

In this section, we study the partial transport problem on Riemannian manifolds where the 
cost is taken to be the square of the Euclidian distance d. Indeed, existence and uniqueness 
of the partial transport has been established by Figalli [2, Remark 2.11]. Therefore, our main 
concern here will be the regularity of the free boundary. In view of the method developed in the 
previous section, we will solely focus on giving sufficient conditions for local semiconvexity of 
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the free boundary (for definitions, etc. regarding optimal transport in the Riemannian setting, 
the reader may e.g. consult [TO]): 

Theorem 4.1. (Semiconvexity) Let M be a smooth n - dimensional Riemannian manifold with 
Riemannian distance d. Consider two non-negative, integrable functions f = fxn, 9 = 9XA 
with £1 n A = and cut(il) Pi A = 0, where cut(Q) is the cut locus of f2. Furthermore, assume 
A is bounded and 4- - convex with respect to Q. Then the free boundary in the partial transport 
problem with cost c := ^ is locally semiconvex. 

Proof. First, by [21 Remark 2.11], the partial transport exists and classical results imply that 
it has the form T m = exp(V x I , m ) for some c - convex function Next, pick a free boundary 
point x £ f2. Then, for e > small, exp^ may be used as a chart between B e (x) and the tangent 
space at x. Since cut(£l) n A = 0, 

cut(B e (x) n (du m n n)) n T m (B e (x) n (du m n n)) = ; 

thus, we may use only one chart (i.e. exp) to project 

(B e (x) n (du m n tt)) u T m (B t (x) n (du m n n)) 

onto the tangent space at x. The cut locus assumption also implies that d is smooth on 0, x A, 
and since £1 n A = 0, we have that it is bounded away from zero. Thus, the level sets of c 



enjoy a uniform ball condition (see Lemma 2.11). Moreover, thanks to the c - convexity of A 



with respect to Q, we may apply Proposition 3.2 to deduce local Lipschitz regularity of the free 



boundary (note that (2.2) follows from Lemma 2.10 which we can apply thanks to our cut locus 



assumption). The Lipschitz regularity combines with the uniform ball property of the level sets 
and readily yields local semiconvexity of the free boundary. □ 

Remark 4.2. We note that one may remove the disjointness assumption in Theorem \4-l\ by 



localizing the problem as in Corollary 3. 7 



If cut(£l) n A ^ 0, then it may happen that T m {x) G cut(x), so the proof above breaks down 
(since the distance function is only smooth away from the cut locus). Nevertheless, there is 
currently some literature available in understanding the proper conditions which ensure that 
this scenario does not happen. Indeed, Loeper and Villani |7, Theorem 7.1] shed some light on 
this issue: given a uniformly regular manifold (see Definition 4.1]) and two densities \i and v 
such that fj, « dvol and v (A) > avol(A) for any Borel set A C M, there is a a > depending 
on \x and v, so that 

inf d(T (x) , cut(x)) > a. 

It is well-known that S n_1 is a uniformly regular manifold (see Example 4.3]), so one may hope 
to utilize Loeper and Villani's theory to develop a full regularity theory for the partial transport 
first on the sphere, and then in a more general setting. We conclude with a family of examples 
which illustrate that neither the "stay away from the cut locus" property nor the c - convexity 
of the target with respect to the entire source are necessary conditions for semiconvexity of the 
free boundary. 

Example 4.1. In what follows, we outline a method for constructing two general densities 
/ = fXUj 9 = 9XA where C S n_1 and A C S n_1 so that A is not ^ - convex with respect 
to fi, yet the free boundary in the partial transport problem is locally semiconvex away from 
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the common region: let f2 = S n and A be a small spherical cap centered around the south pole 
with height, say, X (measured from the south pole). It is not difficult to see that c x (N,A) is 



an annulus (here, N is the north pole); hence, A fails to be c := ^- - convex with respect to $7. 
Assume J" A / is slightly smaller than J" A g (to ensure the existence of a free boundary) , and the 
mass m transported is slightly larger than J A /. Now, enlarge A to a bigger spherical cap with 
height | so that JqI> Ja 5 + e > where e is a small positive constant (this can be accomplished by 
adjusting / and g at the beginning). Then, it can be shown that the spherical cap fii centered 
at the north pole with height ^ (measured from the north pole) is outside the active region: 
if not, suppose fn inU fm > 0, and choose a subset A of fii n U m so that j A f m = 5 < e; 

then choose a subset B of Q \ A so that J B f = 5. In the original partial transport plan, A is 
transported to T m (A), and it is easy to see that the distance between A and T m (A) is bigger than 
2 — | = Now we modify the original plan by replacing the mass in A by the mass in B, and 
from B to T m {A), we can cook up a new transport map thanks to the mass balance condition. 

„ /l_(Z)2 

Next, let 29 be the largest distance on 0, and note 9 < tan(0) = ^ — ^ — = < i. Thus, 

8 

26> < f < ^ < d(A,T m (A)). Therefore, it is not difficult to see the new plan is cheaper than 
the original, contradicting optimality. Thus, the original partial transport problem is equivalent 
to a new one with source 0, \ ill and target A; in the new problem we do not have a cut locus 
issue - this ensures an interior ball condition; moreover, it is not difficult to see that if x is a free 
boundary point away from the common region, c x (x,A) is contained in a cone on the tangent 
space with vertex at x whose opening is strictly less than tt. Thus, we may proceed as in the 



proof of Theorem 4.1 to obtain local semiconvexity of the free boundary away from the common 
region. 
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